Collective Excitations, NMR, and Phase Transitions in Skyrme Crystals 
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At Landau level filling factors near v = 1, quantum Hall fer- 
romagnets form a Slcyrme crystal state with quasi-long-range 
translational and non-coUinear magnetic order. We develop 
an effective low energy theory which explains the presence 
in these systems of magnetic excitations at low energies be- 
low the Larmor gap (A) and which predicts a dramatic en- 
hancement of the nuclear spin relaxation rate by a factor of 
10^. The effective theory predicts a rich set of quantum and 
classical phase transitions. Based in part on accurate time- 
dependent Hartree-Fock calculations of the ordered state col- 
lective excitation spectrum, we discuss aspects of the T—u — A 
Skyrme crystal phase diagram. 
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At Landau level filling factor v — 1/m (with m an 
odd integer), the ground state of a two-dimensional elec- 
tron system (2DES) is an incompressible strong ferro- 
magnet W, i.e., it has complete spin alignment even in 
the limit of vanishing Zeeman coupling strength and it 
has a gap for charged excitations. Quantum Hall ferro- 
magnets have the unusual property IQ-Q] that their topo- 
logically non-trivial spin texture excitations (skyrmions) 
carry charge. This property profoundly affects their 
physics. In particular, the ground state at filling factors 
slightly away from v = 1/m contains a finite density of 
skyrmions and these are expected P-R] to crystallize, at 
least when they are sufficiently dilute. 

In an external field the spin-wave modes of a coUinear 
ferromagnet have an excitation gap equal to the Zee- 
man energy. On the other hand, a non-coUinear magnet 
can have two Goldstone modes [g| and one of these can 
remain gapless in an external field. The Skyrme crys- 
tal state has non-collinear magnetic order ^J. A sin- 
gle skyrmion spin texture has its spins aligned with the 
Zeeman field at infinity, reversed at the center of the 
skyrmion, and has non-zero XY spin components at in- 
termediate distances which have a vortex like configura- 
tion, y,^ The classical (or quantum mean-field) energy 
of a skyrmion is independent of the angle Lp which de- 
fines the global orientation of the XY spin components. 
This extra C/(l) degree of freedom for a single Skyrmion 
leads to broken symmetry in the crystal ground state 
and hence to a spin wave mode which remains gapless in 



the presence of the Zeeman field. The existence of this 
gapless spin mode (presumably in overdamped form in 
a Skyrme liquid state) dramatically alters the low tem- 
perature physics, manifesting itself in both rapid nuclear 
spin relaxation [|| and giant apparent specific heat [Q . 

In this Letter we address the impact of thermal and 
quantum fluctuations on the physics of Skyrme crys- 
tals. We propose a rich zero temperature phase dia- 
gram in which quantum fluctuations destroy the mag- 
netic order of the Skyrme crystal for some values of 
v and Zeeman coupling strength. Where order exists 
in the ground state, we estimate the finite tempera- 
ture Kosterlitz-Thouless phase transition temperatures. 
Our analysis is based on a boson Hubbard model, which 
we argue describes the low-energy physics of this sys- 
tem, and on time-dependent Hartree-Fock approxima- 
tion (TDHFA) calculations of the collective modes of the 
Skyrme crystal which accurately capture the intricate mi- 
croscopic physics of the crystal state and fix the Hubbard 
model parameters. 

The simplest symmetry-allowed orientation-dependent 
coupling between neighboring skyrmions is of the form 
V — JJ2uj) <^os {ipi — fj). For typical Zeeman coupling 
strengths and v not too close to 1, detailed Hartree-Fock 
calculations B show that the Skyrme crystal is a square 
lattice with opposing orientations { ipi = tpj + tt) for 
skyrmions on opposite sublattices. The low-energy effec- 
tive Hamiltonian we propose is motivated by microscopic 
theory from which it follows |T^ that, in addition to the 
positional degrees of freedom, skyrmions carry an integer- 
valued internal quantum number (K) which is the num- 
ber of flipped spins (relative to the maximally polarized 
state of the same charge). Quantum states of definite 
K are related to the classical states of definite orienta- 
tion by [|ll],|2| \K) c>c / dip exp{iK ip)\ip) , so that for each 
skyrmion Ki and (pi are canonically conjugate. Thus our 
effective model is equivalent to a boson Hubbard model 
p3[ in which the boson number on the ith site is mapped 
to the number of fiipped spins in the ith skyrmion. In 
this language the orientation-dependent interaction term 
corresponds to boson hopping and favors long range bo- 
son coherence. 

We limit our attention for the moment to states where 
the skyrmions have crystallized. At low energies, phonon 
excitations of the lattice decouple |I4| from the spins 



and can be treated separately. Because of the Larmor 
gap and the incompressibility gap for particle-hole excita- 
tions, fluctuations in K are the only relevant low-energy 
degrees of freedom, if the skyrmion positions are fixed. 
The energy of an isolated skyrmion in state Ki can be 
written in the form |^ exi = Uxi + g*iJ.BBKi, where g* 
is the host semiconductor g-factor. Non-linearity in the 
function U{K) maps to on-site boson interactions. In 
the Zeeman coupling term on the right hand side of this 
equation, — g = —g*fiBB, plays the role of a chemical po- 
tential for the bosons. (In typical experimental systems 
g ~ 0.015(e^/e^).) Making a quadratic approximation for 
Uk we are led to an effective Hamiltonian describing bo- 
son fluctuations for states with K near (K) = K^{g, v): 



H = C/ ^(i^, - Ky + J^ cos (</?, ~^j). 
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For J ^ U, boson hopping dominates over interactions 
so that the bosons condense into a supersolid phase with 
long range phase coherence. A linear spin wave approx- 
imation for this model yields a gapless Goldstone mode 
with dispersion 



£;(k) = {4UJ [2 - cos{k^a) - coa{kya)]}^^'^ 



(2) 



However, unless 2K'^ is an odd integer, quantum fluctua- 
tions will destroy the classical ordered state and produce 
a gapped 'insulating' state Q when U :$> J. 

The parameters of this model can be fixed using mi- 
croscopic Hartree-Fock calculations Q . From the ground 
state calculation we can evaluate the number of reversed 
spins per charge, (K). In the harmonic approximation, 
the boson interaction parameter U is given by the inverse 
boson compressibility: 
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The boson hopping parameter J can then be extracted 
from the collective mode spectrum which we compute by 
locating poles of the time-dependent Hartree-Fock ap- 
proximation (TDHFA) spin and density response func- 
tions, Xm^ (k -I- G, k -I- G', Lo) . The TDHFA calculation 
is facilitated by simplifications which arise from the re- 
stricted Hilbert space of the lowest Landau level ||l5|,Q . 
The collective mode spectrum in the extended BZ of the 
square lattice Skyrme crystal is illustrated in Fig. m. The 
two gapless modes are the phonon mode and the new 
spin-wave collective mode discussed above. The dashed 
lines show the same modes folded into the smaller BZ 
of the magnetic lattice. As fc — > 0, the higher energy 
of the two-zone folded modes evolves into the Larmor 
mode (global precession of electronic spins around the 
Zeeman axis) which must occur at hujo — g* ^bB. The 
fact that the (new) spin wave mode at k equal to a recip- 
rocal lattice vector of the magnetic lattice is equivalent 



to the k = Larmor mode can be understood in terms of 
the opposite orientations of skyrmions on the two sublat- 
tices. |1J] Because of the inhomogeneous magnetic order 
of the Skyrme crystal state, spin and density response 
functions couple at general points in the BZ and poles 
occur in all response functions at each collective mode. 
Fig. n^ confirms that at low energies the only relevant de- 
grees of freedom are skyrmion orientations and locations, 
and that these give rise to well separated spin-wave and 
phonon modes. As A: — > 0, the phonon mode is lowest in 
energy and has the characteristic fc'^" magnetophonon 
dispersion. We estimate J by requiring Eq. (0) to repro- 
duce the spin-wave velocity of the microscopic TDHFA 
calculations. In Fig. g we plot the boson Hubbard model 
parameters [K^ , U and J) for square lattice Skyrme crys- 
tals at a series of g values and for filling factors v = 1.05, 
ly = 1.1, ly = 1.15, and ly = 1.2. 

Even the T = 0, u — g phase diagram of the Skyrme 
crystal states is very rich. Fig. ^ summarizes conclusions 
(some still qualitative) we have drawn from present calcu- 
lations. The line of primary importance in this phase dia- 
gram is the dark solid line which delimits |l6[ the stability 
region of the square lattice. In our calculations the insta- 
bility of the square lattice is indicated by vanishing long- 
wavelength magnetophonon energies. The square lat- 
tice shear constant is negative for classical electrons and 
this structure is therefore expected [14| to become un- 
stable when skyrmion orientation-dependent interactions 
become weak. The square lattice is stable in the small 
g, large jz^ — 1| region where the skyrmion size is compa- 
rable to the crystal lattice constant. At low skyrmion 
density, the critical value of g is expected to vary as 
|jy _ i|3/2 corresponding iQ to a fixed ratio of skyrmion 
size to Skyrme crystal lattice constant. We have not yet 
completed a systematic survey of possible crystal struc- 
tures outside this region where, in any event, skyrmion- 
disorder interactions will frequently play a larger role 
than orientation-dependent skyrmion-skyrmion interac- 
tions. Inside the square Skyrme crystal portion of the 
diagram our estimates place J/U > 1 except at large 
:/ — 1. We therefore expect to have a ground state with 
translational and magnetic order in the shaded portion 
of the phase diagram which excludes regions where K'^ is 
near an integer, J/U < 1, and we expect magnetic-order 
to be destroyed by quantum fluctuations. This phase 
diagram is subject to quantitative alteration when quan- 
tum fluctuations in the skyrmion positions become large 
at large \iy — 1|. Outside of the square-lattice portion of 
the phase diagram, Skyrmion orientations will be more 
weakly coupled and we expect quasi-long-range magnetic 
order to be rare. 

We can extract a rough upper bound for the Kosterlitz- 
Thouless critical temperature associated with the loss of 
quasi-long-range magnetic order: Tkt = {T^/2kB)ps ^ 
0.008(eVe^) ~ IK. Here p^ = J is the spin stiff- 



ness for the classical XY model obtained when U = 
and the numerical estimate is for typical magnetic 
fields ~ 10 Tesla. This temperature should be com- 
pared with the classical melting temperature of the crys- 
tal which can be estimated using KTHNY |1^ the- 
ory. For elastically isotropic systems with long range 
interactions, Tm = (&^M/47r) where ji is the shear con- 
stant and h the lattice constant of the crystal. Ne- 
glecting the elastic anisotropy of the square lattice 
[Qq,l3, we estimate /i for the Skyrme crystal from 
the long-wavelength dispersion of the magnetophonon 

mode Tnjj{k) Ki {2-k fx(.e^ / e) (ki)^^^. These considera- 
tions lead to melting temperatures a.t g = 0.015(e^/e^) 
of 5.23,3.63 and 2.51 (in units of lO'^ie'^/el) ~ O.IK) 
for 1/ — 1.10, 1.15 and 1.20 respectively. Therefore we 
expect that the melting transition will occur at a lower 
temperature than the magnetic transition. Dislocations 
in the 'tetratic' fluid state frustrate the skyrmion orien- 
tational order, so the spins presumably disorder at the 
same temperature in a single transition. |1J] For classical 
electrons the melting temperature of the electron crys- 
tal is lower than the KTHNY estimate by a factor of 
~ 2; quantum positional fluctuations ||l9|| cause a further 
reduction. Since oricntational order is responsible for 
the positive T = shear modulus of the square Skyrme 
crystal, thermal and quantum fluctuations of the XY or- 
der are likely to lead to an even larger reduction in the 
present case. 

At i^ = 1, spin-relaxation in a quantum Hall ferro- 
magnet is activated because of the Larmor mode gap. 
For |i^ — 1| / 0, however, both phonon and spin-wave 
Goldstone modes of the Skyrme crystal can relax nuclear 
spins. The (spatially averaged) relaxation rate has a Ko- 
rringa temperature dependence. 



and Xph are numerical constants. (For i/ = 1.10 and 
g = 0.010(e^/^) our detailed numerical calculations yield 
Xsw — 13 and Xph — 1-9.) We expect that, in prac- 
tice, disorder pinning of the Skyrmion lattice will gap 
the phonon spectrum and suppress t~^. The linear rela- 
tionship we find between t~^ and |j^ — 1| is in agreement 
with experiment. Using parameter values from Fig. 2, 
Eq. ^implies that for |i^ — 1| ~ 0.1, the spin-relaxation 
rate should typically be ^ IQ^ times faster than for a 
2D Fermi gas system at S = 0, in rough quantitative 
agreement with the experimental results of Barrett et al. 
P,l3. This enormous enhancement in 1/Ti can bring 
nuclei into thermal equilibrium with the electrons and 
dramatically increase the apparent specific heat of the 
2D electron system ^^. We speculate that the sharp 
peak in the apparent specific heat observed by Bayot et 
al. |1^ may be associated with critical slowing down at 
the skyrmion melting temperature which suppresses mo- 
tional narrowing of the nuclear resonance in the well and 
thereby enhances thermal coupling to nuclei outside the 
quantum wells. pQ| If so, the skyrmion melting temper- 
ature at |j^ — 1| ^ 0.2 is ~ 30 — 40 mK, a believable ^ 10 
times smaller than the naive KTHNY estimate. 
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where J7 is the nuclear resonance frequency. X+- has 
contributions proportional |1J] to S{hfl — ej(fc)) where 
ej{k) is one of the collective mode dispersions shown in 
Fig. 1. At long wavelengths neither phonon nor spin- 
wave excitations change the total electronic spin Sz so 
that there is no contribution to T^^ from the G = term 
in Eq. m However both change S^ locally, and contribute 
G 7^ terms to Eq. which can be extracted from our 
TDHFA calculations. We find that ill 
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where t~^ and t")^ are the spin-wave and phonon con- 
tributions to the relaxation rate normalized to the Kor- 
ringa expression for the same 2DES at i3 = 0, and Xgw 
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FIG. 3. T = phase diagram for Skyrme crystal 
states. The value of g at which the square lattice 
shear modulus vanishes has been determined numerically at 
jz^ — Ij = 0.05,0.1,0.15, and 0.2 and interpolated to delimit 
the square lattice stability region. The solid circles indicate 
the variation of Ko = {K) along the square lattice region 
boundary. Non-coUinear magnetic order survives quantum 
fluctuations in the region (approximately) indicated by the 
shading. 
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1. Collective mode energies in e /{e€) units for 



V = 1.10,5 = 0.016(e^/ef), along the direction T - X. The 
shaded region is the magnetic Brillouin zone. At each k a 
mode is labeled by the response function component with the 
largest residue in x(k, k). 



FIG. 2. Parameters of the low energy effective model for 
V = 1.05 (diamonds), u = 1.1 (plusses), v = 1.15 (circles), 
and V — 1.2 (squares) as a function of g/(e^ /el). In the 
right panel the solid symbols give values of the boson interac- 
tion parameter U while the open symbols give values of the 
boson hopping parameter J. The noise in the latter curve 
reflects numerical uncertainty introduced by the fitting pro- 
cedure used to extract the spin-wave velocity from calculated 
collective mode energies. 
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